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Abstract 

In the framework of a QCD relativistic potential model we evaluate the form factors 
describing the exclusive decay B — > ttttIu. The calculation is performed in a phase 
space region far away from the resonances and therefore is complementary to other 
decay mechanisms where the pions are produced by intermediate particles, e.g. in the 
chiral approach. We give an estimate of the contribution of the non resonant channel 
of the order of B(B~ -» ir + ir~fvi) ~ 2.2 x 10" 4 . 



In this letter we shall study the B— meson decays 
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7T + 7T i U£ 
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Vl ■ 
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^From the experimental side these decays are interesting in view of the future programs 
at the 5-factories. For example, some of the preliminary studies on the CP violations at 
these machines [p]] have examined the possibility to extract the angle a of the unitarity 
triangle by the B — > pix non leptonic decay channel. The non-resonant decay mode B — > 3ir 
would be interesting to analyze in this context, as it might provide a significant background 
to the main decay process. While a calculation from first principles is not available at 
the moment, a useful approximation might be the factorization approximation and, within 
this approximation, the decay modes (|1[) and (|2|) would provide the crucial hadronic matrix 
elements needed to compute the relevant amplitudes. In passing we note that there is another 
channel, i.e. the B~ — > 7i°n°£~U£ decay mode, which will not be examined here because it is 
less interesting from an experimental point of view. 

^From a theoretical standpoint semileptonic 5-meson decays with two hadrons in the 
final state represent a formidable challenge as they involve hadronic matrix elements of weak 
currents with three hadrons. They can be studied by pole diagrams, which amounts to a 
simplification because only two hadrons are involved in the hadronic matrix elements. This 
is the approach followed in some papers where these decays have been examined in the 
framework of the chiral perturbation theories for heavy meson decays J|;@- This method 
is based on an effective theory implementing both heavy-quark and chiral symmetry 0,0], 
0, and allows to achieve, for systems comprising both heavy (Q) and light (q) quarks, 
rigorous results in the combined mq — > oo, m q —* limit. However the range of validity of 
this approach is limited by the requirement of soft pion momenta. In the soft pion limit the 
amplitude is dominated by a few tree diagrams with resonances as intermediate states, and 
some clear predictions can be made, but, at least for B decays, the actual phase space is 
relatively large and the phenomenological interest of these predictions is modest. The aim 
of this letter is to examine the decays ([I]) and (|2]) in the framework of a QCD relativistic 
potential model |7j and to extend the kinematical range where theoretical predictions are 
possible. We shall present a detailed analysis of the four form factors relevant to ([!]); for 
reasons of space we shall only give a prediction for the width of the decay (^). We shall not 
include final state interactions in our calculation since no consistent way to compute them is 
presently available. It's clear however that they can modify our numerical calculations ||. 

In two recent papers: || ,|l(| we have presented an analysis of some semileptonic and rare 
B decays into one light hadron employing the relativistic potential model in an approximation 



1 



that renders the calculations simpler. We wish to exploit here this approximation in the study 
of the B — > wlv decays . 



Let us start with a description of the model (for more details see J7J, || and ]TD[). In this 



approach the mesons are described as bound states of constituent quarks and antiquarks 
tied by an instantaneous potential V(r), which has a confining linear behaviour at large 
interquark distances r and a Coulombic behaviour ~ — a s (r)/r at small distances, with a s (r) 
the running strong coupling constant (the Richardson's potential [[[]]] is used to interpolate 



between the two regions^). Due to the nature of the interquark forces, the light quarks are 
relativistic; for this reason one employs for the meson wave function \l/ the Salpeter |12 
equation embodying the relativistic kinematics: 



Vf + ml + yj-Vl + m 2 2 + V{r) 



¥(f) = Aftf(f) , (3) 



where the index 1 refers to the heavy quark and the index 2 to the light antiquark; M is 
the heavy meson mass that is obtained by fitting the various parameters of the model, in 
particular the b-quark mass, that is fitted to the value m& = 4890 MeV, and the light quark 
masses m u ~ = 38 MeV, m s = 115 MeV. The 5-meson wave function ty(f) in its rest 
frame is obtained by solving eq. (^); a useful representation in Fourier momentum space was 
obtained in || and is as follows 

ip(k) = 47T^m B a 3 e~ ak , (4) 

with a = 2.4 GeV -1 and k = \k\ the quark momentum in the B rest frame; this is the first 
approximation introduced in ||. 

The constituent quark picture used in the model is rather crude. There are no propa- 
gating gluons in the instantaneous approximation : the Coulombic interaction is assumed 
to be static. Moreover, the complex structure of the hadronic vacuum is simplified: the 
confinement can be introduced by the linearly rising potential at large distances, but the 
chiral symmetry and the Nambu-Goldstone boson nature of the 7r's cannot be implemented 
by the constituent quark picture. For these reasons, while there are good reasons to believe 
that eq. ([$]) may describe the quark distribution inside the heavy meson, one cannot pretend 
to apply it to light mesons. Therefore pion couplings to the quark degrees of freedom are 
described by effective vertices. 

To evaluate the amplitude for semileptonic decays, it is useful to follow some simple rules, 
similar to the Feynman rules by which the amplitudes are computed in perturbative field 



*Spin terms are not included in V(r), which, for heavy mesons, is justified by the spin symmetry in the 
limit rriQ — > oo. Their neglect cannot be justified for light mesons, which is one of the reasons why one does 
not use the constituent quark picture for the pions. 
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theory. The setting of these rules is the main innovation introduced in |9j as compared to 
0. For the decays ([I]) and (fj) we draw a quark- meson diagram as in fig. |l]and we evaluate 
it according to the following rules: 

1) for a charged pion of momentum p n we write the coupling 

N q N q , 

7 75 , (5) 

where f n = 130 MeV. The normalization factors N q , N q > for the quark coupled to the meson 
are discussed below; 

2) for the heavy meson B in the initial state one introduces the matrix: 



B = 15^*1 , o H^s) — r (6) 

V3 y m q m b + qi- q 2 2m b 2m q 

where m b and m q are the heavy and light quark masses, gf, q% their 4— momenta. The nor- 

f k 

malization factor corresponds to the normalization < B\B >= 2 and / - — ^li^ik)] 2 = 

J (27r) d 

2rriB already embodied in (|6]). One assumes that the 4— momentum is conserved at the ver- 
tex Bqb, i.e. gf + q% = = -B-meson 4— momentum. Therefore q± = (E b , k), q% = (E q , —k) 
and 

E b + E q = m B ; (7) 

3) to take into account the off-shell effects due to the quarks interacting in the meson, one 
introduces running quark mass m b (k), to enforce the condition 



E = ^m 2 (k) + \k\ 2 

for the constituent quarks []; 

4) the condition ml > implies the constraint 



< k < k M ~ ^ , (9) 



on the integration over the loop momentum k 

d 3 k 



(2tt) 3 ' 



(10) 



5) for each quark line with momentum q and not representing a constituent quark one 
introduces the factor 

x G(q 2 ) , (11) 



^By this choice, the average < m&(fc) > does not differ significantly from the value mj, fitted from the 
spectrum, see fl for details. 



where G(q 2 ) is a shape function that modifies the free propagation of the quark of mass m q i 
in the hadronic matter. The shape function 



m G — m q i 
m G ~ <? 2 



(12) 



was adopted in || and [ lOfl , with the the value m G = 3 GeV 2 for the mass parameter; 
6) for the weak hadronic current one puts the factor 

N q N<ff(l - 75) - (13) 



The normalization factor N q is as follows: 




(if q = constituent quark) 



(otherwise) 



(14) 



7) finally the amplitude must contain a colour factor of 3 and a trace over Dirac matrices; 
for the 7T° coupling a further factor ±^= is introduced (the upper sign for coupling to the up 
quark, the lower one for coupling to the down quark). 

This set of rules can now be applied to the evaluation of the hadronic matrix element for 
the decay ([!]), corresponding to the diagram in fig. |l|; the result is : 



B(P) 



J^(q) 



q 1- q 



v P2-q 2 



q 2 



ft (p 2 ) 



Figure 1: The Feynman diagram for B — > tt + tt semileptonic decay. 



r =< 7r + (pi)7r-(p 2 )|^(l - i 5 )b\B-(p) >- 
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lV3 f d 3 k 9[k M - k] ^{k) G[( gi - q) 2 } G[( qi - q - Pl ) 2 } 

X 



4/2 J (2vr) 3 ^E b E q (m q m b + q 1 ■ q 2 ) (?i ~ ?) 2 ~ m \> (<li ~ Q ~ Pi) 2 ~ m l" 
Tr[{4i +m b )(4 2 +m q ) ^ 2 (^i+ 4~ 4\ - m q ») i>x{{- 4x -m q ,)^{\ -75)] . (15) 

The amplitude with a 7r° in the final state: 

4 = < n + (px)7r°( P2 )\ur(l - l 5 )b\B°(p) > (16) 

is obtained from J M (pi, p 2 , p) as follows: 

Jq(Pi, P2, P) = ~ -j= [J^(pi, P2, P) - J^(P2, Pi, p)] ■ (17) 

Following [0 we introduce the various form factors as follows. We put q = p — p\ — p 2 
and we write 

< 7r + ( Pl )n-(p 2 )\uY(l - l 5 )b\B(p) >= 

= i w+ (P1+P2Y + i (pi -p 2 y + irq* + 2h e^ s p a pipp 25 ■ (18) 

It is useful to introduce the following variables: 

s = (pi+p 2 ) 2 
t = {p-pxf 

u = (P-P2) 2 , 

that satisfy 

s + t + u = q 2 + m 2 B + 2m 2 n . (19) 

The form factors h, r, w + are functions of three independent variables. One can 
choose as independent variables s, t ,q 2 or, alternatively, s, Ei, E 2 , where Ei, E 2 are the 
pion energies in the B rest frame. The relations between the two set of invariants are: 

t = m 2 B + m 2 — 2m b Ex 
q 2 = s + m% - 2 m B {E x + E 2 ) . (20) 

Kinematical range is as follows: 

4m 2 , < s < m 2 B 

< q 2 < {m B - v^) 2 



m 2 B + 2ml + q 2 - s A V s ~ 4m ^ < < m| + 2m 2 + g 2 - s A^/g-^4mj 

2 2 ~ ~ 2 + 2 ' ^ 
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where 

A = J(m 2 B -q 2 + s) 2 - 4m| s . (22) 



From dig) one can extract the different form factors by multiplying by appropriate 
momenta. One gets : 

h _9 J » Pa Pl ? P2S (O'W 

s[(t- ml) {q 2 -t) - st] 1 } 

9 . - 2 s p% + 2 (m| - u - t) (gi + p 2 )" J M + (f-u) (gi - p 2 )" J M 



4 s m| - (2m| - u - t) 2 + (u - t) 2 

(Pi J rP2YJ l i 



. (Pl-P2)^J^ t-U 

W- = % — r (25) 

s 2s 



t + u m 2 
2s ~s 



B 



(26) 



The calculation of the trace in (|15|) is straightforward and is similar to those performed 
in and [10] for similar processes. The evaluation of the integral is however more compli- 
cated, because the kinematics is more involved, due to the presence of an extra momentum. 
The integration can be performed numerically, but is time consuming, because, unlike the 
semileptonic decays with one hadron in the final state, where the integration involves one 
variable, here the integration domain is genuinely three-dimensional. The calculation be- 
comes simpler putting the light quark mass m q = in the relevant formulae, which is an 
approximation we perform and is justified by the small value of m q in our model. Similarly 
we put m n = 0. Nevertheless the computation remains huge, since each of the four form 
factors depends on three variables and the number of points needed to have a good accuracy 
is high. 

An important point to be stressed is the kinematical range in which the predictions of 
the present model are reliable. We cannot pretend to extend our analysis to very small pion 



momenta for the following reasons. First of all, as discussed in fllQfl , when \p w \ — >• the 
results of the model become strongly dependent on a numerical input of our calculation, i.e. 
the value of the light quark mass m q . The numerical value of m q cannot be fixed adequately 
because the values of the quark masses were fitted from the heavy meson spectrum, which 
is not very sensitive to m q ( for more details see ). Therefore the value we consider 
in the model m q ~ 38 MeV (or m q = in the present approximation) has considerable 
uncertainties. For large or moderate \p n \ this uncertainty does not affect the numerical 
results: the pion momenta are sufficiently large to render the results insensitive to the actual 
value of m q . For very small \p n \ the numerical results depend strongly on m q , which makes 
them unreliable. This is the first reason to exclude the soft pion limit from the analysis 
also in this paper. A second reason is that, in the soft pion limit, the role of pole diagrams 



(3 



such as those studied in becomes relevant. These diagrams cannot be accounted for by 
the present scheme, which at most can be used to model a continuum of states, according 
to the quark-hadron duality ideas. The low-lying resonances, such as those studied in [0 
should be added separately^. The same should be said about the resonances encountered at 
small s, such as the p- resonance. This resonance is not considered in ||, but is expected 
to play a major role; indeed experimentally one has BTZ(B° — > p~~£ + vt) = 2.5^ o x 10 -4 , 
which shows that this is a relevant piece of the B— decay width into two pions. Therefore 
we assume a lower cutoff s > so, with so = 1 GeV 2 and we expect that the results are not 
affected by the above-mentioned theoretical uncertainties. Since the p— resonance and the 
chiral contributions discussed in are absent in our approach, their contribution should be 
added separately. We expect a large contribution from the p and a tiny contribution from 
the diagrams discussed in since they are significant in a very small region of the phase 
space (see the discussion in 0). 

For s > Sq our model has no similar limitations. By duality we would expect that the 
sum over higher mass resonances can be reproduced fairly well by the continuum model we 
employ here: therefore these higher resonances should not be separately considered to avoid 
double counting problems. It could be observed, in this context, that the failure observed in 



| T0|| at high q 2 (for the B — > 7r semileptonic decay) would correspond, in the present case, to 
the small s, not to the large s region. 

Instead of presenting the form factors as functions of s, q 2 and t we prefer to consider s, 
Ei and E2, the pion energies. In terms of E\, E2 and s the allowed kinematical range is as 
follows (we put m n = 0): 

s < s < m 2 B 

s ^ rriB 

< E 2 < 



2m B ~ 2 

S -±^-E 2 < E x < — (27) 
2m B AE 2 v ' 



In tables 1-4 we present some numerical results for the form factors h(s, Ei, E 2 ), r(s, Ei, E 2 ), 
w + (s, Ei,E 2 ) and W-(s, Ei,E 2 ) in the B~ — > 7r + n~ semileptonic decay. In each table we 
present all the form factors at fixed s and different values of the (E\, E 2 ) pair (s = 1 GeV 2 
in table 1, s = 5 GeV 2 in table 2, s = 10 GeV 2 in table 3, s = 19 GeV 2 in table 4). These 
results should allow to get a quantitative assessment of the numerical relevance of the var- 
ious form factors in the allowed kinematical range. A graphical presentation of the fitted 
numerical output for s = 5 GeV 2 is given in figs. for the four form factors. A different way 

*This is the reason why in |h]] the B* pole of the B — * tt form factor is not reproduced in the \p^\ ^0 
region. 



7 



r form factor 



h form factor 




Figure 2: The relevant form factors for B — ► tt + tt semileptonic decay at s = 5 GeV 2 as a 
function of E\ and i?2- Units are GeV^ 1 for w + in (c), W- in (d) and r in (a), GeV -3 for h 
in (b). 

to present the data is to introduce averaged form factors. We choose to perform an average 
in the pion energies according to the following formula: 

I rm B /2 f {s+m 2 )/(2m B )-E 2 

M = am / „ / „ , £i> 30 > ( 28 ) 

A(s) Js/(2m B ) Js/(4E 2 ) 

valid for all the form factors (/ = h, r, w + , Here A(s) is the allowed area in the 

(-Ei, -E 2 ) plane: 

4 2 

The numerical results we obtain have an average error around 10%. A simple way to present 
the data is by an analytical formula: for example the data can be fitted by the following 
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relation 

f t \ P (s ~ si) {s - s 2 ) . , 

m = (s-s 3 )i + si ' (30) 

a procedure which introduces an average numerical error of ±10%; we stress, however that 
in computing the width we have not used this fit and therefore this further error has not 
been introduced. We also point out that the Breit-Wigner shape of eq. ( |3~0|) is a useful 
parameterization and has no dynamical meaning. 

The values of the coefficients Sk , j3 appearing in (|30|) are reported in table 5 for all the 
form factors of the B~ decay. The form factors are depicted in figs. |3|. 

We observe that due to the limitations of our approach, the kinematical region of validity 
of the present model has no overlap with the soft pion region where pole diagrams, see e.g. 
||, are expected to dominate. Therefore a comparison of our work with the results of these 
models is impossible. 

Let us now evaluate the partial width T(B — > mriu). The relevant formulae to compute the 
width are reported in and we do not reproduce them here. We only give our numerical 
results for the cut-off width (s > lGeV 2 ). Numerically we get 

BK{B- -> 7r + 7r-r^) = 2.2 ( -J^-^ j x 1(T 4 (s > lGeV 2 ) . (31) 

For the other decay channel we have: 



BK(B° -> 7r + 7r°r^) = 3.2 ( n } Vub } n \ x 10" 4 (s > lGeV 

\3.2 x 10 A j 



(32) 



The contribution of the p resonance to these decay modes can be estimated in the present 
model as follows: 

BK{B~ -> 7r + 7r-r^) = 1.2 ( J Vub } n , ) x 10^ 4 . (33) 



p ' V 3.2xl0- 3 



For the other decay channel we have: 



V, 



2 



BTZ{B" - tt+tt^) p = 2.4 l ^f^ l x 10" 4 . (34) 

This latter branching ratio is in agreement with the experimental figure quoted above. 

We can therefore conclude that from an experimental point of view the semileptonic decay 
channel with two non-resonant pions in the final state represents an interesting process with 
a significant branching ratio, of the same order of magnitude of the single pion or the single p 
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(c) (d) 

Figure 3: The averaged form factors for B" — > 7r + 7r~ semileptonic decay. Units are GeV 
for w + in (c), W- in (d) and r in (a), GeV~ 3 for h in (b). 



semileptonic decay mode. It would be nice to find this decay mode in the future experimental 
analysis and to test the present prediction of the QCD relativistic potential model. 
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Table 1: Numerical values of the form factors for several values of Ex , E 2 (in GeV) and 
s = 1 GeV 2 . Units are GeV -1 (r, w + and u>_ ) and GeV -3 (h). 



{Ex, 


E 2 ) 


r 


h 


w + 




(0.14 


2.36) 


- 0.26 


0.53 


5.99 


3.12 


(0.18 


2.07) 


- 0.23 


0.62 


5.89 


2.83 


(0.23 


1.79) 


- 0.16 


0.73 


5.89 


2.44 


(0.28 


1.51) 


- 0.039 


0.88 


5.77 


1.92 


(0.35 


1.23) 


0.18 


1.07 


5.60 


1.28 


(0.43 


0.94) 


0.56 


1.32 


5.41 


0.56 


(0.57 


0.66) 


1.22 


1.60 


5.18 


- 0.019 


(0.85 


0.38) 


2.21 


1.72 


4.20 


0.58 



Table 2: Numerical values of the form factors for several values of E\ , E 2 (in GeV) and 
s = 5 GeV 2 . Units are GeV -1 (r, w + and u>_ ) and GeV" 3 (h). 





E 2 ) 


r 


h 


w + 




(0.54 


2.40) 


2.15 


- 0.76 


- 1.25 


- 13.9 


(0.62 


2.16) 


2.36 


- 0.84 


- 0.31 


- 15.1 


(0.71 


1.92) 


2.60 


- 0.92 


0.85 


- 16.5 


(0.82 


1.68) 


2.88 


- 0.99 


2.27 


- 17.9 


(0.96 


1.44) 


3.22 


- 1.06 


3.98 


- 19.3 


(1.14 


1.20) 


3.62 


- 1.10 


6.01 


- 20.7 


(1.40 


0.95) 


4.08 


- 1.08 


8.29 


- 21.5 


(1.82 


0.71) 


4.52 


- 0.95 


10.5 


- 20.9 



Table 3: Numerical values of the form factors for several values of E\ , E 2 (in GeV) and 
s = 10 GeV 2 . Units are GeV" 1 (r, w + and u>_ ) and GeV" 3 (h). 



(Ei, E 2 ) 


r 


h 


w + 


w_ 


(1.03 , 2.45) 


0.42 


- 0.29 


- 2.71 


- 0.36 


(1.13 , 2.26) 


0.29 


- 0.32 


- 2.71 


- 0.33 


(1.24 , 2.08) 


0.14 


- 0.34 


- 2.68 


- 0.37 


(1.37 , 1.89) 


- 0.044 


- 0.37 


- 2.60 


- 0.49 


(1.52 , 1.70) 


- 0.24 


- 0.40 


- 2.42 


- 0.75 


(1.70 , 1.51) 


- 0.45 


- 0.42 


- 2.11 


- 1.15 


(1.93 , 1.32) 


- 0.63 


- 0.44 


- 1.60 


- 1.78 


(2.23 , 1.13) 


- 0.75 


- 0.44 


- 0.81 


- 2.72 
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Table 4: Numerical values of the form factors for several values of E\ , E2 (in GeV) and 
s — 19 GeV 2 . Units are GeV -1 (r, w + and u;_ ) and GeV -3 (h). 



(Ei, E 2 ) 


r 


h 


w + 


U7_ 


(1.87 , 2.55) 


- 0.81 


0.0041 


- 1.31 


1.73 


(1.94 , 2.45) 


- 0.87 


0.0028 


- 1.37 


1.77 


(2.02 , 2.36) 


- 0.95 


0.0012 


- 1.44 


1.80 


(2.10 , 2.27) 


- 1.02 


0.00092 


- 1.51 


1.82 


(2.19 , 2.17) 


- 1.10 


- 0.0033 


- 1.57 


1.83 


(2.29 , 2.08) 


- 1.18 


- 0.0062 


- 1.64 


1.83 


(2.40 , 1.99) 


- 1.27 


- 0.0096 


- 1.70 


1.81 


(2.51 , 1.89) 


- 1.35 


- 0.014 


- 1.75 


1.77 



Table 5: Numerical values of the parameters appearing in the formula 

f(s) = ^-j 1— JL _ y( s ) i s ari y f the averaged form factors. Units are 

\s — s 3 ) + s 4 

GeV 2 for s\ , s 2 , S3 and S4; f(s) and (3 have the same units. 



form factor 


Sl 


S2 


«3 


s 2 
s 4 


P 


w+ 


- 0.957 


+ 7.15 


+ 3.98 


+ 3.38 


- 1.87 


W- 


+ 1.77 


+ 11.4 


+ 4.20 


+ 3.47 


+ 3.91 


h 


+ 3.25 


+ 19.5 


+ 5.38 


+ 7.93 


+ 0.207 


r 


- 0.446 


+ 9.47 


+ 4.32 


+ 11.5 


- 1.57 
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